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ABSTRACT 

In  all  over  the  world,  land  use  land  cover  changes  ( LULCC )  effect  on  the  water  storage  structures  is  important 
for  proper  management  of  water  resources  and  its  uses.  The  effects  of  LULCC  and  improper  management  systems  have 
played  a  significant  role  for  sediment  transport  movement,  regulation  and  life  expectancy  of  the  reservoir.  For  the 
analysis  of  LULCC  in  Finchaa  watershed  the  geographic  information  system  (GIS)  version  based  SWAT  has  been  used 
to  evaluate  sensitivities  and  patterns  of  LULCC.  To  predict  stream  flow  because  of  LULCC  calibration  and  validation  of 
the  soil  and  water  assessment  tool  were  applied  using  the  compatible  version  of  SWAT  CUP  against  stream  flow  for 
Finchaa  hydropower  reservoir  were  taken  to  estimate  model  performance  on  monthly  basis.  LULCC  involve  all  the 
parameter  which  are  responsible  to  decrease  the  proper  regulation  and  performance  in  ecosystem.  To  predict  stream 
flow  changes  the  analysis  were  performed  on  surface  runoff,  sediment  yield  and  flow  of  ground  water.  It  was  found  that 
the  SWAT  model  predict  stream  flow  due  to  LULCC  reasonably  well  having  R2  values  of  0.71  &  0.74  and  Nash-Sutcliffe 
efficiency  were  found  to  be  0.68  &  0.72  for  calibration  and  validation  of  the  model  resulting  the  best  model  to  predict 
hydrological  process  of  the  basin  or  catchment. 
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1.  INTRODUCTION 

Land  and  water  are  the  two  most  precious  and  fundamental  resources  essentially  required  not  only  for 
nourishment  of  life  but  also  for  the  economic  and  social  development  of  the  country.  The  main  cause  of  LULC 
changes  due  to  hydrological  events  such  as  overland  flow  runoff,  infiltration,  evaporation,  percolation  and  stream 
flow  (Li,  et  al.,  2007,  chakravarti  et  al.,  2014,  2015)  and  also  due  to  human  induced  factors  such  as  deforestation, 
urbanization,  land  grazing,  wood  and  fuel  utilization.  Africa  is  experiencing  rapid  and  substantial  social,  economic, 
climatic  and  environmental  change.  Ethiopia  experiences  persistent  land,  water  and  environmental  degradation  due 
to  localized  and  global  climatic  variances.  The  dynamics  of  land  use  is  one  of  the  major  environmental  problems  in 
Ethiopia  (Berhan,  2010).  Man-  made  reservoirs  usually  satisfy  multiple  objectives  including  flood  control, 
irrigation,  hydropower  generation,  water  supply,  boating,  fishing  and  recreation.  Blue  Nile  Basin  is  the  bigger  area 
in  Ethiopia  having  big  population  density,  land  degradation  and  mostly  required  on  agricultural  economy  (Tsegaye, 
2006).  The  main  factors  of  erosion  and  stream  flow  problems  are  land  use  changes  which  have  been  induced  by 
Finchaa  reservoir.  The  main  purpose  of  this  study  is  to  calibrate  and  validate  the  hydrologic  SWAT  model  based 
on  a  stream  flow  data  of  hydropower  reservoir,  using  Geographic  Information  System  (GIS)  based  version  of  the 
model  to  use  as  efficiency  as  compared  to  other  model  to  access  the  available  resources  such  as  LULCC, 
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deforestation  and  change  of  watershed  management.  With  respect  to  available  water  resources,  any  LULCC  may  affect  the 
hydrological  cycle  and  its  water  balance  terms.  Increase  in  the  intensity  of  rainfall  has  a  significant  impact  on  soil  erosion 
rates  and  increase  the  rates  of  runoff  to  the  reservoir.  The  rain  storm  intensity  cause  significant  erosion  and  sediment 
movement  resulting  the  operation  and  maintenance  cost  increase  and  reducing  life  of  hydraulic  structures  (Tamene,  et  al., 
2005). 

Finchaa  hydropower  reservoir  is  a  highland  area  with  a  severe  soil  erosion  problem  that  drains  to  the  Blue  Nile 
River,  which  may  result  in  reducing  the  capacity  of  water  storage  in  the  reservoir.  The  development  of  Finchaa  dam  and 
Finchaa  Sugar  Estate  caused  the  changes  of  land  use  and  the  rate  of  environmental  degradation  downstream  of  Finchaa 
watershed  in  general.  Therefore,  to  address  the  above  situation,  watershed  management  is  one  of  the  most  important 
approaches,  which  helps  to  reduce  land  degradation,  increase  vegetation  cover,  and  increases  the  productivity  of  the 
watershed  area.  A  suitable  examination  of  the  runoff  yield  of  the  catchment  is  essential  for  management  of  sediment  yield 
and  operation  of  water  resource. 

2.  STUDY  AREA  DESCRIPTION 

The  study  area  is  located  in  Horro  Guduru  Wollega  zone.  East  Wollega,  Oromia  regional  state,  western  Ethiopia 
between  9°46'45"-9°10'30"  North  latitude  and  37o28'30"-37o03'00"  (Figure  1).  Finchaa  sub-basin  is  a  part  of  Blue-Nile 
basin.  The  sub-basin  has  an  area  of  4089  km2.  Finchaa  dam  was  constructed  in  1973  for  the  development  of  economic 
growth  in  the  country  and  it  was  expanded  by  diverting  Amerti  river  flows  into  the  Finchaa  reservoir,  by  constructing  a 
20m  high  earth  and  rock  fill  dam  and  1.57  km  long  diversion  tunnel,  thereby  diverting  Amerti  river  to  reservoir  in  1987.  At 
the  time  the  dam  was  the  largest  hydro-electric  Projet  in  the  country.  The  Finchaa  hydropower  dam  which  has  a  340  meter 
crest  length  and  height  above  lowest  foundation  level  is  20  meter. 

3.  MATERIALS  AND  METHODS 

The  Main  Materials  used  for  Input  Data  Preparation,  Analysis  Were 

•  Arc  SWAT,  SWAT-CUP,  PCP  STAT,  Dew02.exe,  Microsoft  Excel,  DEM,  Meteorological,  Hydrological  map 
and  data 

The  methods  used  here  in  this  research  work  involved  validation,  calibration,  sensitivity  and  uncertainty  analysis. 
SUFI-2  used  as  calibration  and  uncertainty  analysis  algorithms  were  used. 

The  Methodology  of  this  Work  has  the  Following  Components 

•  Collection  of  data,  processing,  input  model,  sensitivity  test,  calibration  &  validation  of  model.  Model  result 
analysis. 

Geographical  Information  System-GIS  based  version  of  soil  and  water  Assessment  Tool-SWAT  were  applied  to 
analyze  overall  methodology.  Finally,  calibration,  validation  and  evaluation  by  appropriate  systems  to  check  the 
model  performance  with  observed  data. 

3.1  SWAT  Model  Description 

SWAT  has  already  been  validated  in  the  regions  of  the  world  for  a  variety  of  applications  in  hydrologic  and  was 
developed  for  the  reason  of  simulation  estimate  the  management  of  land  practices,  sediment  movement  in  large  of  the 
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watersheds  by  changing  soils,  land  use  and  agricultural  conditions  over  extended  time  periods  (Neitsch,  et  al.,  2005). 

SWAT  can  be  used  to  analyze  small  or  large  catchments  by  discrediting  them  into  sub-basins,  which  are  then 
further  sub-divided  for  modeling  purpose.  The  catchment  is  divided  into  a  number  of  sub-basins  which  will  be  divided  into 
hydrological  response  units  (HRUs)  each  having  homogeneous  land  use,  soil  types,  and  management  and  slope 
characteristics.  The  water  balance  in  hydrological  response  units  is  calculated  on  the  bases  of  day  by  day  precipitation, 
evapotranspiration,  overland  flow,  infiltration  and  subsurface  and  groundwater  return  flow. 

SWt  =  SW0  +  I[=1(Rday  -  Qsurf  -  Ea  -  wseep-Qgw)  (3.1) 

Where  SWt=  Final  soil  water  content  in  mm,  SW0=  Initial  water  content  in  mm,  /?day=  Precipitation  on  day  i  in 
mm,  Qsurf =  Surface  runoff  on  day  i  in  mm,  Ea  =  evapotranspiration  on  day  i  in  mm,  wseep=  Water  entering  the  vadose 
zone  from  the  soil  profile  on  day  i  in  mm,  Qgw  =  Return  flow  on  day  i  in  mm,  t  =  Time  in  days. 

3.1.1.  Surface  Runoff 


SWAT  uses  two  methods  for  calculating  surface  runoff  namely,  modified  Soil  Conservation  Service  curve 
number  technique  (USDA-SCS,  1972)and  for  infiltration  method  Green  &  Ampt  were  used  (Green  and  Ampt,  1911).  The 
SCS  techniques  were  applied  in  the  present  investigation,  which  is  a  function  of  permeability  of  soil,  land  use  and  AMC 
(antecedent  moisture  content).  In  SCS  technique,  curve  number  changes  non-linearly  with  soil  moisture  content.  Green  & 
Ampt  technique  requires  sub-daily  precipitation  data  and  estimate  infiltration  as  a  function  of  the  wetting  front  metric 
potential  and  effective  hydraulic  conductivity.  The  SCS  curve  number  equation  is  given  by  the  following  equation: 


(Rj,y  -  la)2 

(Rday  “  la  +  S) 


(3.2) 


Where  Qsurf  =  Rainfall  excess  in  mm  H20,  R  day  =  Depth  of  rainfall  in  day  in  mm  H20, 

la  =  Initial  abstractions,  infiltration  &  interception  in  mm  H20,  S  =  Retention  parameter  in  mm  H20. 

The  retention  variable  changes  spatially  because  of  variation  in  soils  type,  land  cover,  slope  and  soil  water 
content.  The  equations  to  define  the  retention  parameter  are  shown  by  the  following  equation: 


Where:  CN  =  Curve  number  in  days.  Ia,  can  be  written  as  0.2S.  So  that  the  equation  3.2  can  be  converted  into  the 
following  equation: 


Qsurf  = 


(^day  0.2S) 
(Rday+0.8S) 


(3.4) 


The  overland  flow  comes  when  Rday  greater  then  Ia.  The  rate  of  peak  overland  flow  is  the  maximum  overland 
flow  (runoff)  rate  that  occurs  with  a  given  rainfall  event.  The  peak  overland  flow  rate  is  an  indicator  to  predict  sediment 
loss.  SWAT  model  estimate  the  peak  runoff  rate  with  a  modified  rational  method  (Neitsch,  et  al.,  2005)  which  can  be 
express  by  the  following  formulas: 


_  C.i.A 

qpeak  “  17 


(3.5) 
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Where:  qpeak  =  Peak  rate  of  runoff  in  mVs,  C  =  Coefficient  of  runoff,  I  =  Intensity  of  rainfall  in  mm/hr,  A  =  basin 
area  in  km2. 


3.1.2.  Sediment  Component 


SWAT  model  determines  the  amount  of  sediment  component  in  terms  of  the  surface  erosion  and  sediment  yield 
under  every  hydrologic  event  unit  with  modified  universal  soil  loss  equation.  The  quantity  of  stream  sediment  routing  is 
based  on  the  high  amount  of  sediment  which  can  be  transported  from  a  reach  segment,  and  is  a  function  of  peak  channel 
velocity  (Neitsch,  et  ah,  201 1).  The  modified  universal  soil  loss  equation  can  be  express  as  follows: 


Sed  H.8 (^Qsiirf' Qpeak‘ A  hru) 


0.56 


KuSLE  ■  C[JSLE  ’  PlJSLE  ■  LS usLE  ■  ^FRG 


(3.6) 


Where:  Sed  =  Sediment  yield  on  a  given  day  in  tons,  Qsurf  =  Surface  runoff  amount  in  mm  water/ha,  qpeak  = 
Rate  of  peak  in  m3/s,  Ahru  =  Basin  area  in  ha,  KUSLE  =  Erodibility  Value. 


The  deposition  and  degradation  are  the  two  processes  of  sediment  transport  in  the  channel  system  in  the  reach. 
Degradation  observes  when  soil  concentration  is  lesser  than  quantity  of  soil  which  can  be  moved  from  a  reach  segment  and 
deposition  observes  when  soil  concentration  is  high  than  the  maximum  amount.  The  sediment  amount  in  the  reach  can  be 
described  by  the  following  equations: 


sedch  sedch  i  SGd^ep  T  SGd^eg  (3.7) 

Where  sedch  =  Quantity  of  suspended  particles  in  reach  in  metric  tons/day,  sedch,i  =  Quantity  of  suspended 
particles  in  reach  at  the  beginning  of  the  time  in  metric  tons/day,  seddep  =  Quantity  of  deposited  particles  in  the  reach  in 
metric  tons/day  and  seddeg  =  Quantity  of  sediment  re-entrained  in  the  reach  in  metric  tons/day. 

The  Quantity  of  transported  particles  out  of  the  reach  can  be  estimated  by  the  following  equations: 

SGdout  SGdCfr  *  SGdfigp  y  ^  (3.8) 

Where  sedout  =  Quantity  sediment  moved  out  of  the  reach  in  metric  tons/day,  Sedch  =  Quantity  of  suspended 
particles  in  the  reach  in  metric  tons/day,  Vout  =  Quantity  outflow  during  time  step  in  m3/s,  and  Vch=  Quantity  of  water  in  the 
reach  in  m3. 


3.1.3.  Estimation  of  Potential  Evapotranspiration 


The  parameter  of  evapotranspiration  is  estimated  by  Penman-Monteith  equation  which  can  able  to  estimate  all  the 
climatic  parameter  required  by  the  SWAT  model.  The  equation  to  estimate  the  above  parameter  is  given  by  the  following 
equations: 


ET0 


0.408A(fln-G)+yjr^t/2(es-ea) 

A+r(l+0.34(/2) 


(3.9) 


Where  ET0  =  Amount  of  evapotranspiration  Rn  =  Net  radiation  at  the  crop  surface  (  ^  ),  G  =  Soil  heat 


\dayj 


flux  density^— ^ T  =  Mean  daily  air  temperature  at  2  m  height  (°C),  U2  =  Wind  speed  at  2  m  height  j,  es  = 
Saturation  vapor  pressure  in  KPa,  ea  =  Actual  vapor  pressure  KPa,  (es  -  ea)  =  Saturation  vapor  pressure  deficit  KPa,  A  = 


Slope  vapor  pressure  curve 


/  KPa\ 

V  0C  ) 


,  y  =  Psychometrics  constant 


/  KPa\ 

V  0C  ) 
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3.2.  Data  Input  and  Their  Sources 

SWAT  model  is  highly  data  concentrated  model  that  requires  specific  information  of  the  basin.  The  input  data’s 
required  for  this  study  were  Digital  Elevation  model  (DEM),  land  use  and  land  cover  map,  soil  map  and  data,  stream  flow 
data,  and  weather  data.  These  data  were  collected  from  different  sources. 

A  Digital  Elevation  Model  of  30  m  by  30m,  in  the  Grid  format  and  projected  was  used  in  this  study  and  the 
original  DEM  in  geographic  coordinate  system  was  obtained  from  Ethiopian  Ministry  of  water.  Irrigation  and  Energy 
bureau  GIS  Department. 

The  LULC  data  for  this  study  area  was  obtained  from  Ethiopian  Ministry  of  water.  Irrigation  and  Energy  bureau 
GIS  Department. 

The  sediment  input  data  required  by  SWAT  model  to  calculate  stream  flow  includes  the  relevant  hydraulic 
conductivity  properties  such  as  bulk  density,  hydraulic  conductivity  available  water  capacity  of  sediment.  The  parameters 
of  sediment  data  were  obtained  from  Ethiopian  Ministry  of  water.  Irrigation  and  Energy  bureau  GIS  Department. 

Daily  stream  flow  was  obtained  from  Ethiopian  Ministry  of  water.  Irrigation  and  Energy  bureau  hydrology 
Department  and  the  available  weather  data  for  the  selected  stations  were  obtained  from  Ethiopian  National  Meteorological 
Service  Agency.  The  average  monthly  minimum  and  maximum  temperature  patterns  of  different  stations  (1985-2014)  (Fig. 
2  and  3) 

3.3.  Input  Data  Preparation,  Processing  and  Analysis 

After  the  data  was  collected,  an  analysis  of  all  the  collected  data  was  made.  One  of  the  problems  in  hydrology 
especially  in  developing  countries  is  hydrological  data  both  in  quantity  (length  of  record)  and  quality  (standard  of  scientific 
approach).  The  output  of  any  research  depends  highly  on  data  input.  The  acquired  data  were  checked  for  any  outliers  and 
missing  values.  The  missing  meteorological  and  stream  flow  data  were  filled  using  linear  regression  method. 

3.4.  SWAT-CUP 

SWAT-CUP  is  an  interface  that  was  developed  for  SWAT.  SWAT-CUP  is  designed  to  integrate  various 
sensitivity  analysis,  calibration,  validation  and  uncertainty  programs  for  SWAT  using  different  interface.  The  function  of 
interface  is  to  provide  a  link  between  the  input  and  output  of  a  calibration  program  and  the  model.  Using  this  generic 
interface,  any  calibration,  validation,  uncertainty  or  sensitivity  program  can  easily  be  linked  to  SWAT. 

3.5.  Model  Calibration  and  Validation 

The  calibration  is  an  intensive  process  used  to  make  the  suitable  variables  in  modeling  and  an  iterative  process 
that  compares  simulated  and  observed  data  of  interest  (typically  stream  flow  data)  through  parameter  evaluation.  The 
degree  of  accuracy  of  parameter  estimation  was  assessed  by  applying  the  model  to  different  data  set  that  was  not  used  for 
calibration. 

3.6.  Model  Efficiency 

The  performance  of  the  model  to  simulate  the  stream  flow  during  the  calibration  and  validation  periods  has  been 
evaluated,  based  on  the  computed  results  of  the  indicators  and  the  suggested  model  performance  rating  standards.  The 
methods  for  goodness  of  fit  measures  the  model  predictions  were  applied  for  calibration  and  validation  process,  these 
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model  performance  measures  the  variation  in  the  measured  data  that  is  replicated  in  the  simulated  model  results  are 
coefficient  of  regression  (R2)  and  the  Nash-Sutcliffe  efficiency  (NSE).  R2  ranges  from  0.0  to  1.0  with  higher  values 
indicating  better  agreement  (Legates  and  McCabe,  1999)  and  the  value  of  NS  ranges  from  minus  infinity  to  1.0,  with 
higher  values  indicating  better  agreement  (Legates  and  McCabe,  1999). 

R2is  calculated  by  the  following  equation: 

_  [X[L1(qsi-qs)(qoi-qo)]2 
XSUCqsi-qs)2  £”=1(q°i-qo)2 

Where:  qsi  is  simulated  stream  flow  in  m3/s,qoi  is  observed  stream  flow  in  m3/s,  qsis  the  mean  of  the  simulated 
value,  qois  the  mean  of  the  observed  value. 

The  Nash-Sutcliffe  simulation  of  the  model  efficiency  indicates  the  degree  of  fitness  of  the  observed  and 
simulated  plots.  It  is  calculated  as  follows  with  the  same  variables  defined  above: 

i  _  £^i(qoi-qsi)2 
£f=1(qoi-qo)2 

3.7.  Model  Setup 

The  SWAT  model  setup  were  concerned  integrating  the  necessary  input  data’s  with  the  model  such  as  DEM, 
watershed  delineation,  land  cover  and  land  use  map,  sediment  description  and  their  property.  The  basin  delineation 
procedure  contains  5  importent  steps,  DEM,  definition  of  stream,  inlet  outlet,  basin  outlets  selection  and  calculation  of  sub 
basin  variables.  SWAT  model  predict  area  of  hydrologic  cycle  for  every  HRU  and  routes  to  predict  total  loadings  of  the 
catchment.  For  this  study,  20%  land  use  threshold,  20%  soil  and  10%  for  slope  was  used.  After  land  use/  soils/slope 
definition  and  HRU  definition  was  done  a  shape  file  called  ‘Full  HRUs’  were  created.  As  per  the  final  report  from  the 
model  the  basin  was  divided  into  27  sub-basins  which  were  further  divided  into  339  hydrologic  response  units. 

The  shape  file  of  the  prepared  land-use  land  cover  was  given  as  input  for  the  model  data  of  the  SWAT  to  describe 
the  HRU  of  the  watershed.  Consequently,  the  impact  of  each  type  of  LULC  was  considered  in  this  model  to  compute 
stream  flow  in  the  basin.  To  integrate  the  soil  map  with  SWAT  model,  manually  define  was  used.  The  model  uses  slope  for 
the  development  of  HRU  in  addition  to  LULC  and  sediment  input  variables.  Arc  SWAT  allows  the  integration  of  land 
slope  up  to  5  classes  when  defining  hydrologic  response  units.  Table  1  shows  the  slope  classes  of  the  Finchaa  watershed. 

4.  RESULTS  AND  DISCUSSIONS 

4.1.  Sensitivity  Analysis 

It  is  important  to  identify  sensitive  parameters  for  a  model  before  applying  SUFI-2  for  calibration  to  avoid 
problems  known  as  over  parameterization  (Van  Griensven,  et  al.,  2006).  To  know  the  sensitive  variables,  Latin  hypercube 
simulation  one  at  a  time  (LH-OAT)  technique  was  applied  (Van  Griensven,  et  al.,  2006).  Twenty  seven  variables  were 
predict  for  model  parameterization  sensitivity  test,  only  twelve  of  them  were  effective  and  sensitive  parameters  that  are 
responsible  for  monthly  flow  simulation  analysis  (Table  2).  It  was  found  that  the  sensitive  variables  were  responsible  for 
the  model  calibration  and  variable  changes  at  model  iteration  processes.  The  result  obtained  from  sensitivity  test  shows  that 
the  twelve  flow  variables  were  sensitive  to  SWAT  model  which  means  that  the  hydrological  event  of  the  research  area 
mostly  dependent  on  the  action  of  these  variables. 
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4.2.  Stream  Flow  Calibration 

For  the  stream  flow  calibration,  the  comparison  can  be  taken  between  the  observed  and  computed  runoff  for 
thirteen  years  duration  steps  from  1990  to  2002  years  on  the  monthly  basis.  Before  running  for  the  above  time-step  years, 
initially  one  year  of  discharge  data  from  1990  was  used  as  the  warming  time  and  the  rest  of  time  was  taken  for  the  model 
calibration.  The  model  is  calibrated  by  taking  12  variables  which  is  observed  as  most  sensitive  variables  was  taken  for 
measurement  of  discharge  and  coefficient  of  determination  (R2)  and  NES  is  used  to  analyze  the  model  efficiency  for 
observed  and  computed  flows  on  monthly  basis  are  found  to  be  0.7 1  and  0.68  respectively.  The  uncertainties  were  measure 
by  taking  percentage  of  measured  data  bracketed  referred  to  P-factor  by  95  percent  prediction  uncertainty  and  have  been 
calculated  as  0.52.  Similarly  strength  of  the  model  calibration  and  uncertainty  has  been  analyzed  using  the  average 
thickness  of  95PPU  or  R  factor  which  shows  the  average  thickness  of  95PPU  band  divided  by  standard  deviation  of  the 
observed  data  which  have  been  computed  as  0.69. 

The  comparison  between  the  computed  and  observed  flow  for  the  station  during  the  calibration  time  was  shown  in 

figure  4. 

4.3.  Stream  Flow  Validation 

Calibrated  parameters  were  validated  for  the  period  of  nine  years  apart  from  that  one  year  in  was  taken  as  warm 
up  time.  The  validation  shows  the  performance  of  model  for  computed  discharge  in  time  different  from  calibration. 

The  correlation  coefficient  R2  =  0.73  and  the  Nash  Sutcliffe  Efficiency  (NSE)  =  0.72  shows  a  good  agreement 
between  the  observed  and  simulated  values.  The  representation  of  observed  and  computed  stream  flow  for  station  during 
the  validation  period  was  shown  in  figure  5. 

5.  CONCLUSIONS 

The  present  investigation  was  performed  to  predict  stream  flow  change  son  surface  runoff,  sediment  yield  and 
flow  of  ground  water  using  SWAT  model  in  Finchaa  watershed  Ethiopia.  The  evaluation  was  done  in  terms  of  the  impacts 
of  LULCC  on  the  variation  of  LULC  types  and  the  variations  caused  depending  on  the  LULC  types  on  the  major 
components  of  discharge  including  overland  flow  and  vertical  movement  of  water  in  the  ground.  The  physically  based, 
spatially  distributed  and  public  domain  Soil  and  Water  Assessment  Tool  (SWAT)  was  successfully  used  to  simulate  the 
impact  of  LULCCs  on  stream  flow  on  the  hydropower  reservoir. 

The  model  was  carefully  calibrated  and  validated  for  the  Finchaa  watershed  and  represent  better  result  for  the 
performance  evaluation  of  the  model.  Therefore,  SWAT  can  be  utilized  very  well  for  hydrological  simulations  in  the 
selected  catchments  which  have  capabilities  to  analysis  of  hydrological  responses  in  the  watershed.  It  was  observed  that  the 
SWAT  model  predict  the  stream  flow  due  to  land  use  land  cover  changes  reasonably  well  having  R2  values  of  0.71  &  0.74, 
and  NES  were  found  to  be  0.68  &  0.72  for  calibration  and  validation  of  the  model  instead  of  any  other  hydrological  model 
or  algorithm  to  represents  the  hydrological  process  of  the  catchments. 


Table  1:  The  Slope  Classes  of  the  Finchaa  Watershed 


Classes 

Slope  Range 

(%) 

Area 

Land  form 

Ha 

% 

Class  1 

0-8 

Flat  or  almost  Flat 

218468 

76.31 

Class  2 

8-30 

Gentle  Sloping,  Undulating  Plain 

56396 

19.70 
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Table  1:  Contd., 

Class  3 

>30 

Steep  hills.  Very  steep  slopes, 
ridges,  and  mountains 

11436 

3.99 

Table  2:  Result  of  the  sensitivity  Analysis  Parameters  of  Flow  in  Finchaa 
Watershed  (Sensitivity  Parameters  using  SUFI-2) 


Parameters 

NO. 

Name 

Description 

1 

CN2 

SCS  runoff  Curve  number  for  moisture  condition  II 

2 

alpha bf 

Base  flow  alpha  factor  (days) 

3 

GW DELAY 

Ground  water  Delay  (days) 

4 

GWQMN 

Threshold  depth  of  water  in  shallow  aquifer  required  for  return  flow  (mm) 

5 

GW REVAP 

Ground  water  ‘revap’  coefficient 

6 

ESCO 

Soil  evaporation  compensation  factor 

7 

CH N2 

Manning’s  roughness  coefficient  for  main  channel 

8 

CH K2 

Effective  hydraulic  conductivity  of  the  main  channel 

9 

alpha bnk 

Base  flow  alpha  factor  for  bank  storage 

10 

SOL AWC 

Soil  available  water  capacity 

11 

SOL K 

Saturated  Hydraulic  conductivity 

12 

SOL_BD 

Moist  bulk  density 

Figure  1:  Location  of  the  Study  Area 
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Figure  2:  Average  Monthly  Minimum  and  Maximum  Temperature 
Patterns  of  Different  Stations  (1985-2014) 
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Figure  0:  Average  Monthly  Rainfall  Distributions  for 
Different  Stations  (1985-2014) 
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Figure  4:  Calibration  Result  of  Monthly  Observed  and 
Simulated  Flows  of  Finchaa  Watershed 
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FLOW  OUT  12 


Monthly  stream  flosr  Validation  (2003-2011) 
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Figure  5:  Validation  Result  of  Monthly  Observed  and  Simulated  Flows  of  Finchaa  Watershed 
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